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This paper is concerned with general circuits of linear, time-varying, 
positive, two-terminal components. It describes methods of manipulating 
corresponding matrix (vector) differential equations. It uses the manipula- 
tions to derive equations for power and bounds on stability. The bounds 
apply to the exponential factors associated with the basis functions of 
periodically varying circuits. 

I. INTRODUCTION 

Linear time- varying circuits of special kinds have been designed and 
analyzed with notable success. On the other hand, theoretical techniques 
suitable for more general linear time-varying circuits have been develop- 
ing much more slowly. The development of more general techniques can 
be approached in various ways. One can seek to specialize the pure 
mathematics of linear differential equations, in order to discover the 
properties of those equations which can actually correspond to physical 
circuits. Alternatively, one can seek to apply the classical analysis of 
general dynamical systems. As still another alternative, one can seek to 
generalize, for time-varying circuits, concepts, principles, and techniques 
which have long been applied to fixed circuits. 

This paper illustrates the circuit theory approach. After formulating 
matrix (vector) differential equations corresponding to circuits of linear, 
time-varying, two-terminal components, it describes some general meth- 
ods of manipulation. These apply to combinations of time- varying 
matrices and the differentiation operator, and are time-varying counter- 
parts of manipulations applied to constant matrices in the theory of 
fixed circuits. Thereafter, the paper uses the manipulations to derive 
formulas for power, and some bounds on stability. 

The power equations are conventional and reflect the well known fact 
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that time-varying capacitors and inductors can supply energy to a 
circuit. However, they can be derived in a way which illustrates manipu- 
lation of time-varying matrices in particularly simple terms. Further- 
more, the form of the power equations suggests a starting point which 
leads eventually to the bounds on stability. 

It is assumed throughout the paper that all circuit components (fixed 
or varying) are positive. The stability bounds (as derived) assume that 
the circuits vary in a periodic manner. Then the basis functions can be 
arranged as a set of exponentials, each multiplied by a periodically vary- 
ing coefficient (except in singular cases which are the time-varying 
counterparts of fixed circuits whose frequency functions have multiple 
poles). The basis functions are the counterparts of the "natural modes" 
of fixed networks, and they play an equally important role. 

The signs of the real parts of the exponents, in the basis functions, 
determine whether the functions grow indefinitely, or die out. The sta- 
bility bounds derived herein arc upper and lower bounds on the real parts 
of the exponents. The specific bounds depend on the composition of 
the circuit — ■ whether it is composed exclusively of resistors and ca- 
pacitors, resistors and inductors, or capacitors and inductors, or in- 
cludes all three kinds of components. For each composition, there are 
two pairs of bounds, corresponding respectively to the node equation 
and mesh equation (which will be defined in Section II). 

The form of the bounds is illustrated by the following: Let G and C 
be the node matrices of the conductances and capacitances of a period- 
ically varying circuit of resistors and capacitors. Consider the zeros in 
terms of the scalar variable X of the determinant of (G + \C + \C), 
If the capacitances are positive at all times, the matrices can be so 
defined that C is positive definite. Then the zeros of the determinant are 
real, and they vary periodically with time. It will be shown that the 
time averages of the instantaneous maximum and minimum over the 
set of zeros are upper and lower bounds on the real parts of the exponents 
in the basis functions. 

The mesh equations lead to similar bounds in terms of the zeros of the 
determinant of (K — \R + A/?), in which K and R are the mesh mat- 
rices of the stiffnesses and resistances of the capacitors and resistors. It 
is well known, on energy grounds, that a circuit of varying positive 
resistors and fixed positive capacitors cannot be unstable. These bounds 
show that, likewise, a circuit of fixed resistors and varying capacitors 
cannot be unstable. It is true even though the varying capacitors can 
give power gain. 
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Similar bounds can be obtained for circuits of resistors and inductors 
only, by a simple transformation of current and voltage variables. 

The bounds for "RC" and "RL" circuits are at least reminiscent of 
known bounds on the characteristic roots of a dissymmetrical constant 
matrix. Let A be a real matrix. The largest and the smallest of the char- 
acteristic roots of H-'l + A ') are upper and lower bounds on the real 
parts of the characteristic roots of A itself. 

For a circuit of inductors and capacitors only, let S be the node 
matrix of the reciprocals of the inductances, and let C be again the node 
matrix of the capacitances. Consider the zeros of the determinants of 
the two matrices (\C + \C) and (-J5 + \8). Treat the two sets of 
zeros as a single set of numbers. Then the time average of the maximum 
over the set is an upper bound on the real parts of the exponents in the 
basis functions, and the time average of the minimum over the set is a 
lower bound. The mesh analysis leads to similar hounds, except that the 
pertinent matrices are now (\L + XL) and (-\K + \K), where L and 
K are the mesh matrices of inductances and stiffnesses. 

In some ways, the bounds for circuits of inductors and capacitors are 
less satisfactory than those for circuits of resistors and capacitors. When 
the inductors and capacitors are fixed (and positive) the damping is 
necessarily zero. When they are varying the damping may or may not 
be zero. The bounds derived here generally bracket zero damping. Thus 
they do not say whether or not a time-varying circuit of inductors and 
capacitors has any basis functions with damping different from zero. 
They merely say that if the damping is different from zero it is at least 
within the bounds. 

Similar bounds arc easily established for circuits of all three kinds of 
components. However, they tend to be weaker than the bounds for 
circuits of two kinds of components only, in ways which will be explained. 

Some of the derivations have not been completed to the extent of 
proving validity for all singular, as well as normal, situations. What is 
reported here is the result of exploratory studies of time-varying circuits, 
which still leave some details to be filled in. 

The author is indebted to I. W. Sandberg, whom he has consulted 
freely concerning properties of matrices in general and of positive definite 
and semidefinite matrices in particular. The circuit and power equations 
are formulated more carefully by C. A. Desoer and A. Paige. 1 However, 
their analysis tends more to pure mathematics and less to the theory of 
fixed circuits. They do not include the circuit theory type of manipula- 
tions or the stability bounds, which are the primary concerns of this 
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paper. H. E. Meadows has noted somewhat similar, but weaker bounds. 2 
R. A. Rohrer 3 has derived independently the same stability bounds. 
However, he has done so in terms of classical dynamics (generalization 
of the equations of Hamilton and Lagrange) rather than in the circuit 
theory terms used here. 

The material in this paper is organized as follows: All aspects of circuits 
of resistors and capacitors are considered first — formulation of circuit 
equations, power relations, and stability conditions. Then the same 
analysis is shown to apply to circuits of resistors and inductors, by 
changes of variables. Thereafter, circuits including both inductors and 
capacitors are analyzed in a similar way. This appears to be less con- 
fusing than trying to develop the properties of all the kinds of circuits 
simultaneously. 

II. CIRCUITS OF RESISTORS AND CAPACITORS 

2.1 Formulation of Circuit Equations 

Circuit equations can be formulated for linear time-varying com- 
ponents in almost exactly the same way as for fixed components. Basi- 
cally, there are two parts to the formulation. The first defines the be- 
havior of individual components; the second applies rules for combining 
the effects of the various components in a circuit. 

The components with which we are concerned are of the two-terminal 
or one-port type. Fig. 1 represents a typical component, with terminals 
j and A-. Voltages Ej and E k are associated with the terminals. Current 
Iki enters the component through terminal k and leaves it through 
terminal j. (Then /,-* = —Ikj.) 

The role of the component is to perform an operation which inter- 
relates the voltage difference Ek — Ej and the current /*,• . The com- 
ponents with which we are concerned here perform only linear opera- 
tions. We shall be interested sometimes in the operation which transforms 
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Fig. 1 — A two-terminal component. 
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(Ek — Ej) into Ikj and sometimes in the inverse operation. Thus we 
may write 

I ki = \ kj (E k - Ej) 



E k - Ej = A»f '/« 



(1) 



in which A*,- is a linear operation and A*/" is its inverse. 

Sometimes it is more convenient to consider the charge Q k , , which is 
of course related to the current by I kj = Q kj . We shall make extensive 
use of the symbol p to indicate differentiation: 

p=|. (2) 

Then the charge-current relation becomes 

hi = pQlj (3) 

or, inversely 

Qki = f hidt. (4) 

The operations performed by linear resistors and capacitors are dis- 
played in Table I . For each of the two kinds of components, the opera- 
tions are stated in the two inversely related forms, and in terms of both 
/a-/ and Qkj • In a time-varying circuit, the resistance R kj and conductance 
G k j of the resistor and the capacitance C k j and stiffness K k j of the ca- 
pacitor may be functions of time. 

When the coefficients vary with time they must be written in proper 

Table I — Linear Operations Performed by Resistors and 

Capacitors 

Notation: E k — Ej = voltage across component 

Ikj = current through component 

Qkj — charge delivered to component 

hi = pQkj 

Resistors: G ki = conduct auce, R kj = resistance 

Rk, = Okf 1 

hi = G ki (E k - Ej), Qkj = f G ki {E k - Ej)dt 
(Ek - Ej) = Rkjhs , (Ek - Ej) = Rt/pQki 

('a pari tars: Ckj = capacitance, K k j = stiffness 

Ckj = AV 
Qki = C k j(Ek - Ej), h-j = p<'kj(E k - Ej) 

E k - Ej = KkjQkj , E k - Ej = A'*, f hjdt 
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relationship to the differentiation operation p. For example 

pC kj (E k - Ej) = i [C kj (E k - Ej)} * C kj p(E k - Ej). (5) 

Throughout the paper we shall be critically concerned with the "non- 
commutability of p and time-varying coefficients." 

A circuit is formed by interconnecting a number of components, for 
example as illustrated in Fig. 2(a). The interconnections may be repre- 
sented by the corresponding linear graph, as in Fig. 2(b). The interac- 
tions between the various components are determined by Kirchoff 's two 
laws. The voltage law says that the same voltage E k can be assigned to 
each node (graph vertex) fc, in forming the voltage differences for all 
components connected to k. Then the sum of the voltage differences 
around any mesh (graph cycle) must be zero. The current law states 
that the sum of the currents into any node must be zero. These remarks 
are exactly the same whether or not the components vary with time. 

We shall consider separately the two different forms of circuit equa- 
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Fig. 2 — A typical circuit: (a) circuit diagram; (b) linear graph. 
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tions which are commonly used — the node equations and the mesh 
equations. For the node equations one node is chosen as datum, and the 
excitation (forcing function) is described as currents fed into other nodes 
as in Fig. 3(a). Then the node voltages (relative to the datum) are 



(a) 




NODE = DATUM (l No = - 2 I Nk , E = o) 




Fig. 3 — Node and mesh currents and voltages: (a) node analysis; (b) mesh 
analysis. 

related to the currents by a vector differential equation. The procedure 
is exactly the same for time-varying circuits as for fixed circuits, pro- 
vided one is careful to preserve the correct order of the differentiation 
operator p and time-varying coefficients. 
The vector equation may be written 

I N = (G + P C)E N . (6) 
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Here E N and I N are column matrices, or vectors, whose elements are 
the voltages and excitation currents associated with the various nodes. 
If only certain of the nodes are externally accessible, for excitation cur- 
rents, the elements of I corresponding to the other ("internal") nodes 
are simply constrained to be zero. 

G and C are square matrices defining the specific relation between I N 
and E N . Their elements are easily determined from the instantaneous 
conductances and capacitances of the resistors and capacitors of the 
circuit, and the usual relations of elementary circuit theory apply even 
though the conductances and capacitances vary with time. 

It follows that G and C, at each instant of time, have the properties 
of the matrices usually associated with fixed networks. Thus G and C 
are symmetrical, and if the conductances and capacitances of the com- 
ponents are nonnegative, G and C are positive definite or semide finite. 
The symmetry and the nonnegative character of the matrices, together, 
lead to an important part of the specialization in circuit theory, rela- 
tive to the usual pure mathematics of differential equations. 

When the differentiation operator p appears in front of a matrix, it 
signifies differentiation of each element of the matrix. When it is fol- 
lowed by a matrix product, it signifies differentiation of each element in 
the single matrix equal to the product of matrices. Thus (6) is merely 
a compact way of writing 

In = GE N + — Q N 

Q N = CE n = (<zy) (7) 

d n (d 
dt QN= \dt qi 

For the mesh equations, meshes are chosen in a somewhat arbitrary 
way, and the excitation is described as voltage generators inserted in 
the meshes as in Fig. 3(6). The meshes correspond to cycles in the linear 
graph, and the number of meshes is the maximum number of independ- 
ent cycles permitted by the topology of the graph. The mesh currents 
are related to the excitation voltages by a vector differential equation. 
As before, the procedure is the same for time-varying as for fixed cir- 
cuits, provided p and time-varying coefficients are written in the proper 
order. 

For a network of resistors and capacitors, it is more convenient to 
use charges in place of currents. Then 
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E M = (K + Rp)Q M 

(8) 
Im = pQm ■ 

Here E M , Q\t , Iu are column matrices or vectors. E M and I M are not 
the same as E N and I N , although they are related to them in a quite 
complicated way. 

R and K are square matrices defining the specific relation between 
E M and Q M . Their elements are determined from the instantaneous 
resistances and stiffnesses of the resistors and capacitors of the circuit, 
and the usual relations of elementary circuit theory again apply. R and 
K are related to G and C of (6), but in a quite complicated way. The 
elements of R and K are not simply reciprocals of elements of G and C 
nor are matrices R and K the inverses of and C. 

Like and C, R and K have properties of matrices usually associated 
with fixed networks. They are symmetrical, and if the resistances and 
stiffnesses of the components are nonnegative, the matrices are positive 
definite or semidefinite. 

A further point should be mentioned. If (8) is to be valid for all 
circuits, the constants of integration implicit in the Q kj , as defined in 
(4), must be consistent in the following sense: They must be such that 
the indefinite integral in (4) can be replaced by a definite integral, say 



f hjdl, 



with t the same for all kj. This is the same as requiring that all Qkj must 
be zero at some one time t . 

When the condition is not met, the superposition theorem can be 
invoked to express the complete circuit relations in two parts. The first 
assumes that all capacitors are completely discharged at time t , and 
relates the charges accumulated at t > t to the generator (mesh) voltages 
at t > to. The second starts with the actual charges at t = t and deter- 
mines their later values in the absence of generator voltages. When the 
capacitances vary with time, the initial charges may have a much more 
important effect than is possible in a fixed circuit. This will be discussed 
in physical terms in Section 2.4.1. 

2.2 An Algebra of p and Matrix Coefficients 

This section introduces manipulations of a sort which we shall use 
extensively in this and later papers. In the manipulations it is convenient 



2584 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 19G3 

to use both p and a dot over a symbol to indicate differentiation with 
respect to time : 

px = X = -3- X. (9) 

Generally (but not quite always) the dot will be used for rates of change 
of coefficients and p for differentiation of primary variables (such as 
voltages or currents) or products of coefficients and primary variables. 
By way of introduction, consider the following scalar expression: 

pax = -7- (ax) = ax + npx = (a + ap)x. (10) 

Suppose £ is a principal variable and a is a time-varying coefficient. 
Then a may be regarded as a linear operator which multiplies x by a 
function of time. In the same way, p is an operator which differentiates 
x, pa is an operator which multiplies x by a and differentiates the prod- 
uct, and ap is an operator which differentiates x and multiplies the 
derivative by a. 

Equation (10) may be said to state a commutation rule, which may 
may also be stated as an "operator identity" (both sides of which are 
operators) : 

pa = d + ap. (11) 

Thus p and a commute without change if and only if a is constant, so 
that a = 0. 

The concept is easily extended to more complicated combinations, for 
example 

(a + p)(b + p) = ah + b + (a + b) p + p 2 . (12) 

An algebra of this sort, in p and scalar coefficients, is useful for the 
manipulation of scalar differential equations. It is a principal tool in, 
for example, Ref. 4. 

For present purposes we need to extend the concept to an algebra of 
p and matrix coefficients, as a tool for manipulating vector differential 
equations. Suppose X is a matrix variable and A is a matrix coefficient. 
As a first example, it is easily established that 

pAX = AX 4- ApX = (A 4- Ap)X. (13) 

The relation follows at once from the fact that each element of AX is a 
sum of terms, each of which is a product of one element from A and one 
from X. Equation (11) can be applied term by term and the results can 
then be sorted into contributions to AX and ApX. 
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Corresponding to (13) is the 'operator identity" 

pA = A + Ap. (14) 

The concept is quickly extended to more complicated operations. Some 
"operators equations" are collected in Table II, together with some 
familiar algebraic matrix identities. In the table 

A ' = transpose of A 

A' 1 = inverse of A (15) 

U = unit or identity matrix. 

Constant scalar and matrix coefficients commute with p without 
change. Time-varying coefficients do not. { But of course matrix factors, 
constant or not, generally do not commute with each other.) This may 
be regarded as the most important distinction between the theories of 
linear time-varying and fixed circuits. If it were not for the difference 
in the commutation rules, most of the familiar techniques applied to 
fixed circuits would apply directly to time-varying circuits. As it is in 
fact, the more complicated commutation rules lead to numerous com- 
plications, as we shall see. 

2.2.1 Matrices of Order One 

Much of this paper is concerned with scalar quantities (for example, 
net input power) derived from vector circuit equations. For some pur- 

Table II — Some Matrix Relations 

A, B, X, Y = matrices, <p = a scalar 

In relations which involve inverses, pertinent matrices are assumed to be square 
and nonsingular. 

SOME ALGEBRAIC IDENTITIES 

X(A + B)Y = XAY + XBY 
ipX\ = Xipll = XYip 

(AT)' = Y'X', (XY)- 1 = F-yr-j 

SOME OPERATOR IDENTITIES 

p A = A + Ap, p*A = A + 2Ap + Ap 1 

pAp = Ap + Ap' = — pA 4- p 2 A 
pAip = >i>A + <ppA, pAe* = e?(<pA + pA) 

DERIVATIVES OF SOME MATRIX FUNCTIONS 

pX* = XX + XX * 1XX 

pXYZ = XYZ + XYZ + XYZ 

£p(A'-') = -xx-\ pX-* = -x-ixx-* * -{x-^yx 
P x' = (pxy 
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poses, scalars may be represented by matrices of order one, and the 
scalars in question will be derived in that form. Certain operations will 
be used repeatedly in this connection. 

Suppose W and V are n X 1 column matrices and Y is an n X n 
square matrix. Then 

YV = a column matrix 

(16) 
W YV = a matrix of order one. 

A matrix of order one is always symmetrical, for there are no off -diago- 
nal terms to interchange in forming the transpose. Thus, using a trans- 
pose rule from Table II, 

W'YV = (W'YV)' = V'Y'W. (17) 

Certain special cases are particularly important. 

If Y 1 = Y and W, V are column matrices, 

W'YV = V'YW. (18) 

The differentiation rules in Table II require 

V (\"YV) = V'YV + V'YV + V'YV. (19) 

Applying (18) gives 

if y' = Y and V is a column matrix 

V'YV = V'YV = § v(V'YV) - \ V'YV. (20) 

finally, suppose the transpose of }' is the negative of Y. Then (17) 
requires 

if Y l = — Y and V is a column matrix 

V'YV = V'Y'V = - V'YV = 0. (21) 

(If a quantity is equal to its negative it must be zero.) 

2.3 Power and Stability in Terms of the Node Equation 

2.3.1 Instantaneous and Average Powers 

We return now to the circuit equation of the nodal analysis: 

Ts = (G + V C)E N . (22) 

The power P supplied to the circuit by the excitation currents is 
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P = Y,E»J»k = E N 'I N . (23) 

k 

Multiplying the circuit equation by E', 

P = E N ' (G + V C)E N . (24) 

Expanding in terms of identities in Table II gives 

P = E N 'CE N + EsCE N + E N 'CE N . (25) 

Applying (20) to the last term leaves 

P = Es'(G + \C)Es + \v(E N 'CEs). (26) 

The power relation confirms, in circuit equation terms, what must be 
expected on physical arguments. Thus \p(E N l CE N ) is the rate at which 
energy is being stored electrically in the capacitors. Then Etf l GE N is the 
rate at which energy is being dissipated in the resistors, and hE N l CE N 
is the rate at which energy is being removed from the circuit by what- 
ever means are used to vary the capacitances. (Recall that increasing a 
capacitance decreases the stored energy per unit charge.) 

The average power P is frequently of interest as well as the instan- 
taneous power P. For the average over a finite interval, say h to h , 
integration of (20) gives 

Ave P = Ave \E N '(G + \C)E N \ + %^^>J[' • < 27 ) 

When E N 'CE N is bounded at all times, the last term approaches zero 
as (f-2 — h) approaches infinity. Thus, for long time averages, 

if E n 'CEn is bounded at all times 

P = Ave [E N '(G + \C)E N ] . (28) 

2.3.2 Linear Transformations on E,\ and I N 

The power equation suggests rearranging the current equation so as 
to emphasize the matrix (G + \C). Operator identities in Table II 
yield 

/,, = [((? +^C) + h( V C + C V )\E N . (29) 

It is now time to introduce a linear transformation of a sort which 
will be used extensively in this and later papers. In particular, let 

E N = NE 

(30) 
/ = N'ls . 
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Linear transformations of this sort are of course standard means for 
diagonalizing matrices. They are also well known as means for generat- 
ing equivalent circuits, when components are fixed. 5 Their appropriate- 
ness for time-varying circuits is by no means obvious, and they do in 
fact lead to serious complications not encountered in connection with 
fixed circuits. As usual, the complications stem from the commutation 
rules. 

Multiplying (29) by N l and using (30) gives 

t = [N'(G + \C)N + \{ V N l CN + N'CNp) 

+ h(N'CN - N'CN)}£. 

Note that 

(N'CN - N'CNY = -(N'CN - N'CN). (32) 

When C is positive definite, there is a transformation N such that: 

N'CN= U 

(33) 
N\G + hC)N = -D. 

Here U is the unit matrix of suitable order and D is a diagonal matrix. 
Defining D with negative sign simplifies the later discussion. The exist- 
ence of a suitable transformation matrix, at each instant, follows from 
elementary circuit and matrix theory (for positive components). When 
C is only positive semidefinite, (6) can be transformed into a new equa- 
tion, in fewer dimensions, with a positive definite C. One such procedure 
is outlined in the Appendix. 

With the negative sign in the second equation of (33), the (diagonal) 
elements d k k of D are identical with the zeros X fc of the determinant of 
matrix (G + \C + XC). 

dkk = Xfc , 

(34) 

det (G + \C + \ k C) = 0. 

In our applications, because G and C are functions of time, the X*. are 
functions of time. If the circuit components are always nonnegative, 
the Xfc are all real. 

When the transformation N is fixed by (33), (31) and (32) become 

/ = (-D + pU + J)tS 

J = ^N-'N - (N-'N)'} (35) 
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The power equation (26) becomes 

P = E'l = ft'N'l = ft'I = -E'Dft + &(&$)• (30) 

(The product ft'. J ft = because J 1 = — J and ft is a column matrix.) The 
power equation could have been transformed directly, but we shall find 
it informative to have also the transformed current equation (35). 

2.3.3 Some Bounds on the Basis Functions 

Carrying the analysis only a little further establishes some interesting 
bounds on the basis functions of circuits which vary periodically. This 
subsection outlines a derivation, but simplifies the argument by means 
of some somewhat restrictive assumptions. The next subsection reviews 
the derivation and removes most of the restrictions. 

The basis functions are counterparts, for time-varying circuits, of the 
familiar natural modes of fixed circuits. In node terms, they are a set of 
linearly independent solutions of (22) for Es with I N = 0. Thus, if E a is 
a vector basis function, 

- (G + pC)E„. (37) 

When I N = 0, P = and (26) becomes 

= E a \G + \C)E a + ±p(E a 'CE a ). (38) 

Also, when I N = 0, 7 = and the transformed equations (35) and 
(36) become 

= (-D+ Up + J)ft a 

= ft a 'Dft„ + \ V (Ejft a ) (39) 

E„ = Nft. . 

If the circuit has n degrees of freedom there are n basis functions in 
the set. They may be chosen in many different ways, but each choice is 
a linear transformation on every other choice. 

If the circuit varies periodically, G and C in (37) and (38) vary peri- 
odically [and also D and J in (39)]. It is well known that the basis 
functions of a linear differential equation with periodically varying co- 
efficients can usually be so chosen that they are exponentials with 
periodic coefficients.* Any exceptions are singular cases which we can 
best take care of in the next subsection. 



* See the discussion of the Floquet-Poincare 1 theorem in a text on differential 
equations, such as pages 78-81 in Coddington and Levinson. 6 
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Thus for a periodically varying network, we can use 

E a = #„exp (s a t). (40) 

Here s a is a constant, the exponential is a scalar factor, and II a is a 
periodically varying vector. For the purposes of this subsection, we can 
best restrict ourselves to circuits for which s„ and H a are real. The re- 
striction will be removed in the next subsection. 

We shall find it convenient to replace H a by another periodic vector 
F, , related to II a by 

H a = Fj (41 ) 

in which 6 is an arbitrary periodic real function of I (with the same 
period as the time- varying circuit components). In terms of F„ , 

E 9 = F. exp {sj. + 0). (42) 

Because the exponential is simply a scalar factor, (42) and the last 
equation of (39) require 

% = f a exp (s,* + 6) , x 

(43) 
F a = NF a . 

For any matrix A , 

pAe ( "° l+e) = [exp (s a t + *)][(«, + 6) A + pA). (44) 

Using this relation in (39) and then cancelling out the exponential 
factor gives 

= [(«, + 6)U - D + J]F„ + pK 

(45) 
= F a '[{s a + 6)U - D]F„ + hp(F a 'F a ). 

Now S a does not remain bounded over an infinite time interval. Hence 
averaging the second equation of (39) does not eliminate the second 
term. On the other hand, F a is bounded at all times, and thus the second 
equation of (45) implies: 

- Ave F a '[(s a + 6)U - D\F e . (46^ 

Since f\ , D and 6 are all periodic, while s„ is constant, the long time aver- 
age is the same as the average over any one period. 

The matrix (sa + 6)U — D is diagonal. That is, all its elements are 
zero except on the main diagonal, where the typical element is 

om = s„ + — X* 

(47) 
Det (0 + \C 4- \ k C) = 0. 
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Then multiplying out the matrix product in (46) gives 

= Ave £ (*, + 6 - \ k )F ak 2 (48) 

k 

in which F ak is the Ath element in the column matrix F a . 

When all the quantities are real (as assumed), condition (48) cannot 
be true if all the coefficients (s a + 6 — X*) are positive at all times, or 
if all are negative at all times. We can use this circumstance to set bounds 
on the exponent s a of the basis function. The arbitrary function 6 has 
been introduced as a means of strengthening these bounds, as explained 
below. 

Let us assume temporarily that the X* are at all times distinct. Fig. 
4(a) illustrates a plot of the X* over one period of their periodic varia- 
tions. Let \ u be the largest X* , and let d M be a choice of such that 

X M — &u = constant. (49) 

Fig. 4(b) illustrates a plot of the corresponding X* — 8 M . Every coef- 
ficient (s„ -f- 6 — Xt) in (48) will now be positive unless s„ is no greater 
than \ M — 6 M . 

Now 6 M is actually determined uniquely by (49). Because M is 
periodic, the average of 6 M over any period is zero, and hence also the 
long-time average. Then, when \ M — 6 M is constant 

Ave (\m — 6 M ) = X., r — 6 M — Ave \ M 



•ji = X M — Ave X, 



(50) 



Thus all terms in (48) will have positive coefficients (s, + 6 M — \ k ) 
unless 

.% ^ Ave \„ . (51) 

Exactly the same sort of procedure leads also to a lower bound on s, . 
Let X„, be the smallest \ k , and 6„, a choice of 6 which makes X m — 6 m a 
constant. Fig. 4(c) illustrates the corresponding X* — 0„, . The bound on 
s„ can be written 

s a ^ Ave X m . (52) 

2.3.4 Recapitulation, Discussion, and Removal of Restrictions 

This subsection reviews the derivation of the stability bounds, and 
removes most of the restrictions. 

Our diagonalization of the matrices G + \C and C assumed that C is 
positive definite. The Appendix describes how an equation with a posi- 
tive semidefmite C can be transformed into one of lower dimensionality 
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Fig. 4 — Characteristic roots: (a) X; (b) X - 0, with to make maximum = con- 
stant; (c) X — 0, withfl to make minimum = constant. 
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with a positive definite C. More exactly, the transformation is such that 
the components of the transformed current and voltage vectors may be 
divided into two autonomous parts. If n and m are the order and rank 
of the original matrix C, in components of the transformed current and 
voltage vectors are related by a differential equation of our standard 
form, in m dimensions, with a positive definite C matrix, and m linearly 
independent basis function subject to our stability bounds. The remain- 
ing n — in components of the current and voltage vectors are related by a 
purely algebraic equation, which raises no questions of stability. 

In our derivation of the bounds on s* , we restricted ourselves to 
periodically varying circuits such that all the basis functions (corre- 
sponding to a particular choice) are exponentials with periodic coeffi- 
cients. In singular instances such a choice is not possible. However, for 
periodically varying circuits, the basis functions can all be exponentials 
with coefficients which are at most polynomials in t with periodic coef- 
ficients. Furthermore, the coefficients are polynomials in t only when 
more than one basis function has the same coefficient Sk in the exponential 
factor exp (s*£). Out of the set of basis functions corresponding to a 
single Sk , one can always be assigned a periodic coefficient //* . This is 
sufficient to establish our bound on St , without regard for the possible 
existence of other basis function with the same s* . 

We also assumed that the constant s* and the periodic function ll k 
are real. The assumption is not actually necessary. With real circuit 
components, complex basis functions can be chosen in conjugate pairs. 
Then equally valid choices are the real and imaginary parts of the com- 
plex functions. Thus we can write, for a complex exponential basis 
function, 

(X. + iY a )e ( " + *" u = Z la exp (s a t) + iZ 2a exp (s a t) 

Zi a = X a cos uj — Y„ sin w a L (53) 

Zia = X, sin u> a t + Y a cos u a t. 

We can now use either Z\ a exp (s a t) or Z 2a exp (s a l) as a basis function, 
in place of H„ exp (s„t). The only difference is that the Z„'s are not pe- 
riodic, as is H a . However, the Z„'s are bounded at all times, and so are 
F a and /•'„ , which are now defined by 

Zi„ or Z>„ = F„e 

(54) 
/''„ = NF a . 
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Furthermore, the boundedness of F a is all that is required, for passing 
from (45) to (46), provided the average in (46) is a long-time average. 
Thus the real part s a , of a complex coefficient s a + iu, , also obeys our 

bounds. 

An alternate proof uses the complex basis function itself, replaces E„ 
by its conjugate in (38) and (39) and sorts out real and imaginary parts 
of the equations. 

The exponents corresponding to circuits of constant resistors and 
capacitors are necessarily real. However, when both the resistors and 
capacitors vary with time, the exponents may be complex. A number of 
specific examples are known, including examples cited by Desoer and 
Paige 1 and by Meadows. 

Finally, we assumed that the zeros \ k of the determinant of (G + 
\C + XC) were distinct at all times. This is not necessary, provided we 
make a simple change in the statement of our bounds on s„ . Suppose 
the various X fc 's crisscross, as in Fig. 5 (which may be contrasted with 
Fig. 4a). The variables X A , and X m are now defined as the instantaneous 
maxima and minima over the set of variables Xi , • • • , X„ . Otherwise 
the bounds are the same as before. 

Thus, for general circuits of periodically varying positive resistors and 
capacitors, if s a is the real part of the coefficient in the exponential factor 
associated with a basis function, 




Fig. 5 — X jW and X m when A*'s crisscross. 
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Ave A,„ g 8, ^ Ave X.„ (55) 

A,„ , X.w = instantaneous min. and max. 

over the set X x , • • • , X„ . 

Det [G + \C + X fc C] = 0. 

When the circuit components are constants, the X fc are constant, and also 
each .% is exactly equal to one of the X* . Why do the averages of the X* 
only furnish bounds when the components vary with time? 

The answer lies in the joint implications of the vector and scalar 
equations 

= [(*+ 6)U -D + J\F a + V F a 

- Ave{A'[(*+ 0)U - D]F a \ 

J = N~ l N - (N-'NY 

./' = -J 

[equations (45) and (46) of Section 2.3.3]. It will he sufficient to con- 
sider only the second-order case, for which (56) represents the following 
collection of equations 

= (* + & - XiA + K + JnPa 

- (,% + 6 - \, )/v, + F a -< - J n F aX (57) 

= Ave [(s, + 6 - X,)/V + (s, + 6 - X 2 )F ff2 2 ]. 

(J u — J»-> = and ./ji = — J 12 because J' = —J.) 

With constant coefficients, J v < = 0, /''„ is constant, and one suitable 
solution is 

(s a + — Xi) = 0, to satisfy the 1st eq. 

f\., = o, to satisfy the 2nd eq. (58) 

(«„ + 8 — Xi) and F a2 = 0, to satisfy the 3rd eq. 

With time-varying components, making (s„ + 6 — Xi) zero in (57) 
leaves 

= L + JJ*A 

= (Xi - X, )/''„, + #„ 2 - ./i-F ffl (59) 

= Ave [(X, - X, )/Vj. 
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These equations are (at least usually) incompatible when J n ^ 0. Thus 
the nondiagonal matrix J, which appears only when components are 
time-varying, is what destroys the identity of the exponents s„ and the 
determinant zeros \ k (as here defined). Recall that J stems from the 
commutation rules applied to combinations of E l , p, C, and E. 

When the resistors are time-varying but the capacitors are fixed, C = 
and (G + |C + XC) becomes (G + \C), which has been studied 
extensively in connection with fixed networks. When C is positive 
definite and G is positive definite or semidefinite, none of the zeros X* 
can be positive, and the same is true of their averages when they are time- 
varying. Then our bounds exclude any positive s a in the exponential 
factors exp (s a t) associated with our basis functions. But an unstable 
basis function (defined as one which grows indefinitely) requires a posi- 
tive s a . Thus our bounds confirm the stability of circuits in which only 
the resistors vary. 

2.3.5 Comparison with a Known Property of Constant Matrices 

The bounds on s„ are at least reminiscent of known bounds on the 
characteristic roots of a dissymmetrical constant matrix. Consider the 
roots s„ of 

Det (A + s„U) = (60) 

in which A is a dissymmetrical matrix. For the closest parallel to our 
analysis assume that the characteristic roots are all real. The equation 
may be rewritten as follows: 

Det (S + J + s a U) = 
S = h(A + A'), S l = S (61) 

./ = \{A - A 1 ), J' = -J. 
When the determinant is zero, there are nonzero vectors X such that: 
(S + / + s a U)X = 0. (62) 

Because J 1 = —J, this implies the scalar equation 

X'(S + s a U)X = 0. (63) 

Diagonalizing S leaves 

X (s a - \ k )X k 2 = 

(64) 

Det (5+ X*tf) = 0. 
The J term in (62) excludes (in general) an X in which all elements 
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are zero except one. Then s„ ^ X* , but the largest and smallest \ k are 
upper and lower bounds on .% . Simple changes in the analysis establish 
the same bounds for the real parts of complex characteristic roots. 

2.4 Stability in Terms of the Mesh Equation 

Section 2.3 dealt exclusively with the node equation. However, the 
mesh equation can be manipulated in almost exactly the same way. 
Recall the mesh equation [(8) of Section 2.1], 

E u = (K + Rv)Qm . (65) 

Note the order of operations Rp here, and contrast it with pC in the 
node equation. 

Define a scalar M by 

M = X) QmhEmh = QmE m . (66) 

A: 

While M does not have the dimensions of power, it has many of the 
mathematical properties of the power function P which we associated 
with the node equations. Multiplying E M by Q M ' and applying operator 
identities gives 

M = Q U \K - \R)Qu + MQm'RQu). (67) 

Compare this with the power equation (26) in Section 2.3. The appear- 
ance of —\li in M, as opposed to -\-%C in P, reflects the difference of the 
order of operations in Rp and pC. The quantity \(Qm'RQm) does not 
represent stored energy, but it is mathematically similar to the stored 
energy function \E N 'CE n . Both these quantities may be regarded as 
Lyaponov functions. 

Proceeding from here exactly as in the analysis of the node equations 
leads eventually to similar, but not identical, bounds. The basis func- 
tions are now mesh charges Q a instead of node voltages E a . When circuits 
vary periodically they can be chosen as exponentials with periodic 
coefficients : 

Q a = H„exp (s a t). (68) 

The bounds on the s a (or the real parts when complex) may be written 

Ave X m ^ s„ ^ Ave X« 

X m , ^m — instantaneous min. and max. 

(69) 
over the set Xi , • • • , X„ 

Det [K - \R + \ k R] - 0. 
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Now the coefficients s a must be identical in the exponentials associated 
with the node and mesh analyses. This is because the ratios of charges 
and voltages remain bounded. Thus the two sets of bounds apply to the 
same set of s. . Then one can form a single pair of bounds by choosing 
the lesser of the two upper bounds and the greater of the two lower 
bounds. 

When the capacitors are time varying but the resistors are fixed, R = 
and (K - \R + \ k R) becomes the familiar (K + \ k R) of the theory of 
fixed circuits. Then the bounds require that such a circuit cannot be 
unstable (with positive circuit components) even though the time-vary- 
ing capacitors give power gain. 

2.4.1 A Complication in Degenerate Special Cases 

In the discussion of the mesh formulation we have ignored a degener- 
ate special situation, which complicates a more nearly general analysis. 
The complications arise when there is a node within the network to which 
two or more capacitors are connected, but no resistors. 

Consider first the simplest example, as illustrated in Fig. 6(a). Here, 
two capacitors are connected in series between nodes j and A', and no 
other components are connected to their common node c. If the capaci- 
tances are constant, one can simply replace the two capacitors by a single 
equivalent. When they are time varying, the substitution may have to 
be more complicated. 

The reason is briefly as follows: Suppose there is a positive charge on 



V—*lk — g — #— <! 

/ \ +Q, -Qa \\ 

(a) 

v^ — © — fl 

/ \ \\ 

(b) 

Fig. 6 — Illustration of a degenerate special case: (a) time-varying capacitors 
with a common node to which no resistors are connected; (b) the Thevenin equiva- 
lent. 



LINEAR TIME-VARYING CIRCUITS 2599 

one capacitor and a negative charge on the other. Then there will be a 
ratio of charges (at any one time) which will produce a zero voltage dif- 
ference between terminals j and k. If the ratio of capacitances varies 
with time, the required ratio of charges varies with time. Conversely, the 
only combination of constant charges which yields zero voltage difference 
(Efc — Ej) at all times is zero charge on each capacitor. Furthermore, if 
there is initially a positive charge on one and a negative charge on the 
other, both charges cannot be reduced to zero through nodes j and k 
alone. On the other hand, when the capacitances vary periodically there 
is a ratio of constant charges such that the voltage difference varies 
periodically and with zero average. 

For the mesh analysis one can invoke Thevenin's theorem, and replace 
the two capacitors by an equivalent capacitor and a periodic, zero aver- 
age voltage source, as illustrated in .Fig. 6(b). When there are several 
capacitors and no resistors to a single node, or a combination of resistor- 
less nodes, several Thevenin voltage generators are called for. A general 
characteristic of circuits which can lead to this sort of complication is a 
resistance matrix R which is singular (positive semidefinite instead of 
positive definite). 

It follows that the initial charges may be much more important in 
time-varying than in fixed circuits. For example, suppose there is only 
one resistorless node and only two capacitors connected thereto. In a 
fixed circuit, the effects of initial charges will eventually die out except 
for the voltage at the single resistorless node. In a time-varying circuit, 
however, the varying Thevenin voltage may produce voltages at all 
nodes forever. 

For the node analysis, one can ignore these complications by simply 
not eliminating the resistorless nodes. When there are resistorless nodes, 
the conductance matrix G is positive semidefinite. When G has order n 
and rank m (and C is positive definite), n — m of the basis functions have 
zero damping. This can be proved by applying to G, with only minor 
alterations, the analysis applied to a semidefinite C as outlined in the 
Appendix. The Thevenin voltages of the mesh analysis set up undamped 
oscillations within the circuit (s a = 0) and correspond to the undamped 
modes of the node analysis. The validity of the derivation of our stability 
bounds, which does not preclude a singular G matrix, implies that the 
bounds will include zero when G is positive semidefinite. 

One way to avoid these complications is to avoid ideal capacitors. If 
a resistor is connected across each capacitor, to represent the leakage 
through any actual component, there are no resistorless nodes. 
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IN. CIRCUITS OF RESISTORS AND INDUCTORS 

The theory of circuits of resistors and capacitors can be applied to 
circuits of resistors and inductors by interchanging currents and voltages 
and also the mesh and node formulations. Since currents and voltages 
are interchanged, one needs variables which are related to voltages in 
the same way that charges are related to currents. Corresponding to node 
voltage Ek , let 

E k = «i (70) 

or, inversely 

** = f E k dt. (71) 

If an inductor is connected between nodes j and k, E k — Ej is the 
voltage across the inductor, and $ fc - $; is proportional to the flux 
linkages within it. The fact that the definition of** in (70) or (71) leaves 
undetermined a constant of integration reflects the fact that constant 
flux linkages produce no voltage across an inductor. 

Table III states the linear operations performed by inductors, in two 
inversely related forms, and in terms of both Ek — Ej and ** — <£>, . 
The inductance L kj and its reciprocal Ski can vary with time, provided 
the order of the differentiation symbol p and the coefficients is preserved. 

For resistors and inductors, the mesh equation corresponds to the 
nude equation (6) of our previous circuits, and is 

E u = (R + P L)Im. (72) 

The mesh voltage and current vectors Eu and I M are defined as before, 
and also the matrix R of the mesh resistances. Then L is the mesh matrix 
of the inductances of the inductors. If one properly chooses constants of 
integration stemming from (71), the node equation is 

Table III — Linear Operations Performed by Inductors 



I ki = current through component 
Ek — Ej = voltage across component 
** — */ = variable related to flux 
E k = p$k 

Lkj = inductance, Skj = reciprocal inductance 

Ski = Lkj' 1 
# t - * y = Lkjhj , Ek - Ej = pLtjhj 

hi = S*/(** - */), ht = Skj S (Ek - Ej)dl 
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Is = (S + Gp)$„ 

(73) 

E N = 1&N 

and corresponds to the old mesh equation (8). Here vector In represents 
the node currents as l)efore, and G is again the matrix of node conduct- 
ances. The elements of vector <$> N are the <$* and S is the node matrix of 
the reciprocals of the inductances. 

The formulation of (73) requires an assumption regarding constants 
of integration, implicit in the definition of the $* , exactly like the as- 
sumptions regarding the Qkj in the formulation of (8). 

One can now proceed exactly as before to obtain power equations 
and bounds on the exponential factors associated with the basis func- 
tions of periodically varying circuits. 

Section 2.4.1 described degeneracies which occur in the mesh analysis 
of circuits of resistors and capacitors which have resistorless nodes. The 
counterparts for circuits of resistors and inductors occur in the node 
analysis of circuits with resistorless meshes. 

IV. CIRCUITS WHICH CONTAIN BOTH CAPACITORS AND INDUCTORS 

4.1 Circuits of Capacitors and. Inductors Only 

For circuits which contain both capacitors and inductors, but no re- 
sistors, the node equation contains the capacitor and inductor terms 
from our previous equations (6) and (73). Thus 

In = S$ N + pCE N 

(74) 
E N = p$N 

or, replacing E N by p$jt , 

In = (S + pCp)* N . (75) 

As before, S is the node matrix of the reciprocals of the inductances, and 
C is the node matrix of the capacitances. 

To obtain the input power P from the excitation, multiply by E N ', 
which is the same as <j,v' 

P = *AS + pCp)** . (76) 

Manipulating in terms of operator identities, again identifying $> N with 
Eh , and using the symmetry of matrices S and C gives 
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p = -$*/&* + \E N 'CE N + $?(**'**) + MEnCEj,). (77) 

The first two terms are the rates at which energy is removed from the 
system by the means used to vary the inductances and capacitances. 
The last two terms are the rates of increase of the electrical energy stored 
in the inductors and capacitors. 

The basis functions <J» ff are solutions of 

= (S + pCp)*,. (78) 

The corresponding power equation is 

= -pir'fo + \E a l CE a + hp&S*.) + hp(E a 'CE a ) 

E. = p*. . 

If the circuit components vary periodically, we can define a bounded 
function F„ by 

*, = F„exp (8<t+ 6). (80) 

There will be a corresponding function F„' defined by 

E a = F.'exp (s a t+ d). (81) 

Because E a = p$ a , FJ is related to F. by 

FJ = (s, + 6)F a + F a . (82) 

Under reasonable circuit conditions (which exclude, for example, dis- 
continuous changes in inductances), F„' is bounded, as well as F a . Using 
(80) and (81) in (79), eliminating the exponentials and averaging gives 

= Ave F a '[(s, + 0)* ~ h>W* + Ave F a 't[(s. + 6)C + \C\FJ. (83) 

If both S and C are positive definite, all matrices can be diagonalized. 
This merely requires a different transformation on F„ and F a ' . Thus let 

F. = Nja , 1'V - NJlV (84) 

and choose N a and NJ in such a way that 

N.'SN. = U, N a "CN a ' = U 

(85) 
N.'SN. = 2D., N„"CN a ' = -2D/. 

Here £/ is again the unit matrix and D„ and D a ' are diagonal matrices. 
After transformation, (83) can be written 

= Ave E {(«* + * - *»)A*' + («• + * - ^-)^V 2 }. (86) 
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Here, F ak and F ak ' are elements of the column matrices /''„ and F a ', and 
A* and X*-' are the corresponding (like-rowed) diagonal elements of D a 
and DJ . As before, all quantities can be made real (provided /''„*• , P„k 
are only required to be bounded, not necessarily periodic). Then s a can- 
not be so large that all terms in the sum are positive, or so small that- 
all terms in the sum are negative. 

If one makes no use of the implicit relation between /'„ and F a ', one 
can consider the two sets of constants A* and A/ as simply two parts of a 
single set. Then bounds on s a can be obtained exactly as before. Since 
the A* depend only on C and the A*' only on S, while C and S represent 
different circuit components, the two sets of bounds are likely to criss- 
cross as in Fig. 7. Thus, including also the relation of A*, and A*' to S and 
C [defined by (85)], the bounds may be written 

Ave (min. over A) ^ .% ^ Ave (max. over A) 

set A = Ajk's and A/'s 

(87) 
Det (-££+ \ k S) = 

Det (|C + \ k 'C) = 0. 

The mesh analysis differs only in the specific quantities involved. The 
mesh circuit equation is 

E M = (K + pLp)Q u 

(88) 
Im = pQti ■ 

As before, K is the mesh matrix of stiffnesses and L is the mesh matrix 
of inductances. 

The input power from the excitation is 

P = -IQu&Q* + \IuLIm + MQmKQm) + hv(I«LIu). (89) 

The bounds on s a may be written 

Ave (min. over A) ^ s a ^ Ave (max. over A) 

set A = Aa-'s and At"s 

(90) 
Det (-hK + Aa/v) = 

Det (£L+ \kL)= 0. 

4.2 Circuits of all Three Kinds of Components 

For circuits which contain capacitors, inductors and resistors, the 
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Fig. 7 — Characteristic roots of (- * S + \S) and (§ C + XC). 

node equation may be written 

Is = (S + Gp + pCp)<S> N 

E N = p*Ar • 
The corresponding power equation may be written 
p = -frjfoj, + \E N l GE N + \E N l CE N 

+ ip(tX) + ±p(E»'CE N ) 



(91) 



[92) 
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in which the new term, \E N 'GEs , represents the power dissipated in the 
resistors. 

The stability hounds may be written 

Ave (min. over X) ^ s a ^ Ave (max. over X) 

set X = Xa.'s and X*"s 

(93) 
Det (-£5+ \ k S) = 

Det (G + ±C + X*'C) = 0. 

Note that Xa-' is the same as X*' for a circuit of capacitors and resistors 
only [see (55)]. On the other hand, the present Xa- do not depend at all 
on the resistors. 

The corresponding mesh equation is 

E M = (K + Ry + pLp)Q„ 

Iu = pQm ■ 

Then the power equation may be written 

P = -IQuKQu + lIu'RI* + \ImLIm 

+ \p{QmKQ M ) + MImLIm). 

Finally, the stability conditions may be written 

Ave (min. over X) ^ s„ ^ Ave (max. over X) 

set A = Xa-'s and X/'s 

(96) 
Det {-\K + \ k K) = 

Det (R + \L + Xa-'L) = 0. 

4.3 Discussion 

Let us return temporarily to circuits of capacitors and inductors only. 
Suppose all the components are constant (and positive). Then it is well 
known that the basis functions (natural modes) are undamped. Since 
»S and C = 0, etc., when the components are constant, all Xa- and Xa-' in 
(88) or (90) are zero. Then the upper and lower bounds must be zero. 
Thus, in the special case of fixed components, our bounds collapse onto 
the actual s a . 

Suppose the inductors are fixed but the capacitors are time-varying. 
Then, for example, in the bounds (88), Xi = but generally Xa-' ^ 0. 
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The zero coefficients X fc will affect the maximum or minimum of the com- 
bined set X at any times when all A*' have the same sign. 

When components vary periodically, C and S, etc., vary periodically 
and C and S, etc., have zero averages. Then any nonzero elements in 
C and S, etc., must be positive some of the time and negative some of 
the time. This suggests (but does not prove) that the bounds on s a in- 
clude zero. It is confirmed by some quite different analysis, which is 
outlined in a short paper in this issue of the B.S.T.J. 7 The paper points 
out that the s„'s of nondissipative circuits occur in pairs of the form 
+s a , — s„ . Thus when any s a has a nonzero real part, another s a will 
have the negative of that real part, and any true bounds are consistent 
therewith. It also follows that the bounds described above can be 
tightened (for nondissipative networks) by using the bound closer to 
zero and its negative. 

The basis functions of a circuit of capacitors and inductors may all 
have zero damping even though the components are periodically time- 
varying. This remark is supported by well known properties of the 
Mathieu-Hill equation, which corresponds to the one-dimensional special 
case of our vector equations. A weakness of our bounds is that they do 
not give meaningful sufficient conditions for basis functions with (posi- 
tive or negative) damping different from zero. The bounds will (at least 
almost) always be different from zero, but by bracketing zero they will 
not exclude it. 

Suppose now the circuit includes resistors as well as inductors and 
capacitors. Consider first the bounds (93) derived from the node equa- 
tion. If the node conductances (represented by matrix G) are sufficiently 
large, every characteristic root A*' will be more negative than every 
characteristic root X fc (which does not involve G). Then the lower (more 
negative) bound on s a will be set by the V and the upper bound by the 
X fc . If the conductances are further increased, the lower bound will be- 
come more negative, but the upper bound will remain unchanged. On the 
other hand, increasing the conductances will usually increase the damp- 
ing of all the basis functions. 

The same remarks apply to the bounds derived from the mesh equa- 
tion except that mesh resistances replace node conductances. Thus when 
resistors add substantial damping in circuits containing both capacitors 
and inductors, our upper (more positive) bounds may be too weak to 
have much significance. 

Our node analysis assumes that both S and C are positive definite, 
although the initial formulation may make one or both positive semi- 
definite. As before, however, the procedure can easily be modified for 
positive semidefinite matrices. Recall that two different transformations 
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are used to diagonalize 8 and (', in the derivation of our hounds. When 
both are positive semidefinite, they can be separately transformed to 
positive definite submatrices bordered by zeros. Because constant transfor- 
mations can be used (in accordance with the Appendix) , the transformed 
S and C retain the derivative relationship to the transformed S and C. 
Then the number of characteristic roots X* , X*' depends on the ranks of 
S and C. When only capacitors are connected to some nodes, the mesh 
analysis will involve some zero damped Thevenin voltages; when only 
inductors appear in some meshes, the node analysis will involve some 
Thevenin currents. 



AIM'KNMX 



Positive Semidefinite Capacitance, Matrices 

When C is only positive semidefinite, one. can transform equation (6) 
as follows: The rank m of C is now less than the order n. There exist 
transformation matrices Ni such that 



HtfNi = 




c 



(07) 



The zeros represent submatrices of zero elements, bordering the com- 
plete matrix with n — m rows and n — m columns of zeros. Then C is 
an in X m positive definite submatrix. 

These remarks apply to time-varying as well as to fixed circuits. 
Furthermore, the matrix Ni can be so chosen that it is constant, provided 
the time-varying capacitances are always positive ( >0, not ^0). In 
fact, Ni can be chosen as a matrix of 0's and l's only. A detailed demon- 
stration is beyond the scope of this paper; briefly it derives from the fact 
that the rank of C is determined by the topology of the capacitor part 
of the circuit, without regard to the values of the capacitances (provided 
they remain >0, so that the topology cannot be a function of time). 

Because iV\ = 0, Ni'p = pNi, and (6) can be transformed (as for 
fixed circuits) into 



E = Njbi , h = Nil. 



(08) 



A further transformation, A^2 , is now easily found, of the form: 



U 

Nn 1 V 


ft 


U N v , 

u 


= 


Go 

G 



(90) 
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Here U represents a unit, or identity submatrix of appropriate order, 
and N12 is an m X (n — m) submatrix. Because of the zero submatrix 
in the upper right-hand corner of the complete matrix, 



U 
N 12 l 





c 


U N12 



= 




c 



(100) 



Thus the capacitance matrix is unaffected by the second transformation. 
The elements of the N l2 portion of JV 2 may be time-varying. However, 
because the time-varying elements combine with no nonzero elements of 
G, AT 2 'pC = pNz'C. Thus the second transformation yields 



/ = 



D? SMS *0*- 



We can now divide the components of the transformed current and 
voltage vectors into two parts, as follows: 

/„ = (!aE„ , in n — m dimensions 

/ = (0 -\- pC)E, in m dimensions. 



The first equation is algebraic, and its diagonalization is routine. The 
second equation is in our standard form (6), with C positive definite, 
and it may be transformed further in accordance with (29) to (36). 

We arrived at this formulation by starting with a one-to-one cor- 
respondence between modes and components of the current and voltage 
vectors, and then transformed the corresponding differential equation. 
Desoer and Paige 1 arrive at the same conclusion by starting with more 
general choices of the current and voltage components, taking account 
of the circuit topology. 
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